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Abstract 

We consider a communication scenario where a source communicates with a destination over a directed layered relay network. 
Each relay performs analog network coding where it scales and forwards the signals received at its input. In this scenario, we 
address the question: What portion of the maximum end-to-end achievable rate can be maintained if only a fraction of relay nodes 
available at each layer are used? 

We consider, in particular, the Gaussian diamond network (layered network with a single layer of relay nodes) and a class 
of symmetric layered networks. For these networks we show that each relay layer increases the additive gap between the optimal 
analog network coding performance with and without network simplification (using k instead of N relays in each layer, k < N) 
by no more than log(iV/fc) 2 bits and the corresponding multiplicative gap by no more than a factor of (N/k) 2 , asymptotically (in 
source power). To the best of our knowledge, this work offers the first characterization of the performance of network simplification 
in general layered amplify-and-forward relay networks. Further, unlike most of the current approximation results that attempt to 
bound optimal rates either within an additive gap or a multiplicative gap, our results suggest a new rate approximation scheme 
that allows for the simultaneous computation of additive and multiplicative gaps. 

I. Introduction 

Analog network coding (ANC) extends to multihop wireless networks the idea of linear network coding (TJ, where an inter- 
mediate node sends out a linear combination of its incoming packets. In a wireless network, signals transmitted simultaneously 
by multiple sources add in the air. Thus, each node receives at its input a noisy sum of these signals, i.e. a linear combination 
of the received signals and noise. A communication scheme wherein each relay node merely amplifies and forwards this noisy 
sum is referred to as analog network coding |2j, ||3). 

The rates achievable with ANC in layered relay networks is analyzed in (3), BJ. In (3), the achievable rate is computed 
under two assumptions: (A) each relay node scales the received signal to the maximum level possible subject to its transmit 
power constraint, (B) the nodes in all L layers operate in the high-SNR regime, where miring; PR,k > V<M — 1, ■ ■ ■ ,L for 
5 > and Pr^ is the received signal power at the fc th node. It is shown that the rate achieved under these two assumptions 
asymptotically (in source power) approaches network capacity. The authors in [4] extend this result to the scenarios where the 
nodes in at most one layer do not satisfy these assumptions and show that achievable rates in such scenarios still approach the 
network capacity as the source power increases. 

However, requiring each relay node to amplify its received signal to the upper bound of its transmit power constraint results 
in suboptimal performance of analog network coding in general, as we show in |6)-|[8). Further, even in low-SNR regimes 
amplify-and-forward relaying can be capacity-achieving relay strategy in some communication scenarios, [5|. Therefore, it is 
desirable to analyze the performance of analog network coding in general layered networks, without the above two assumptions 
on input signal scaling factors and received SNRs. 

Analyzing the performance of analog network coding without such assumptions however, results in a computationally 
intractable problem in general HI, (6). Therefore, to compute the maximum achievable ANC rate in general communication 
scenarios, in |7| we establish a result that significantly reduces the computational complexity of this problem. Further in |8), 
we propose a greedy scheme to bound from below the optimal rate achievable with analog network coding in general layered 
networks. For general layered relay networks, these two results allow approximating the optimal ANC rates within a constant 
gap from the cut-set upper bound asymptotically in the source power. These results not only allow us to exactly characterize 
the optimal ANC rate in a class of symmetric layered networks that cannot be addressed using existing approaches under the 
assumptions A and B, but allow us to do so in a computationally efficient manner. 

In this paper we introduce another approach to reduce the computational complexity of approximating the maximum 
achievable ANC rate in general layered networks. The proposed approach is based on the notion of network simplification, 
introduced in |9] to characterize fraction of the capacity of the Gaussian TV-relay diamond network when only k out of 
N available relay nodes are used. Previously, in cooperative communication literature as in pO[-fT2[, the notion of network 
simplification has been used in the amplify-and-forward relay networks to characterize the achievable rates. However, such prior 
work used network simplification in a restricted sense (selecting the best single relay node among N relays) and considered 
simple network topologies (the source communicating with the destination via a single relay node). In contrast, we provide the 




Fig. 1. Layered network with 3 relay layers between the source V and destination 't'. Each layer contains two relay nodes. 

optimal ANC rate characterization in much general communication scenarios where any number k of relay nodes in a layer 
are used (k > 1) and the source communicates with the destination over any number L of layers of relay nodes (L > 1). In 
this sense, to the best of our knowledge, this is the first work to characterize the performance of network simplification in 
such general layered amplify-and-forward relay networks. 

We show that in the Gaussian iV-relay diamond network (layered network with a single layer of relay nodes L — 1) and a 
class of symmetric layered networks, the network simplification allows us to maintain achievable rates within small additive and 
multiplicative gaps of the maximum ANC rate achievable when all relays in a layer are used. Along with our previous results 
in 0' ®' where we establish that analog network coding is a capacity achieving strategy in some communication scenarios, 
the results in this paper establish that in those scenarios using fewer nodes in each layer, we can still maintain achievable rates 
within small additive and multiplicative gaps to the capacity. Further our results indicate that in general layered networks, the 
maximum achievable ANC rate can be tightly approximated with much lower computational complexity than with existing 
schemes. 

Organization: In Section [IT] we introduce a general wireless layered relay network model and formulate the problem of 
maximum rate achievable with ANC in such a network. Section III addresses the performance of network simplification in 
the Gaussian iV-relay diamond network and computes additive and multiplicative gaps between the maximum ANC rates 



achievable when N and k(k < N), relays are used. In Section IV we consider a class of symmetric layered networks and 
compute additive and multiplicative gaps between the optimal ANC rates obtained when all relays in every layer are used and 
when only a fraction of such relays are used. Section [V] concludes the paper. 

II. System Model 

Consider a (L + 2)-layer wireless network with directed links. The source s is at the layer '0', the destination t is at the 
layer 'L + V, and the relay nodes from the set R are arranged in L layers between them. The l lh layer contains ni relay nodes, 



, j " — M. An instance of such a network is given in Figure [l] Each node is assumed to have a single antenna and operate 
in full-duplex mode. 

At instant n, the channel output at node i,i G RU {t}, is 



E 



hjiXj[n\ + Zi[n] 



-oo < n < oo, 



(1) 



where Xj[n 



is the channel input of the node j in the neighbor set Af(i) of node i. In ([TJ, hji is a real number representing 
the channel gain along the link from node j to node i. It is assumed to be fixed (for example, as in a single realization of 
a fading process) and known throughout the network. The source symbols x s [n],—oo < n < oo, are i.i.d. Gaussian random 
variables with zero mean and variance P s that satisfy an average source power constraint, x s [n] ~ Af(0, P s ). Further, {^[n]} 
is a sequence (in n) of i.i.d. Gaussian random variables with Zi[n] ~ Af(0, a 2 ). We also assume that are independent of the 
input signal and of each other. We assume that the i th relay's transmit power is constrained as: 



E[xt[n]]<Pi, 



-oo < n < oo 



(2) 



In analog network coding each relay node amplifies and forwards the noisy signal sum received at its input. More precisely, 
a relay node i at instant n + 1 transmits the scaled version of yi[n], its input at time instant n, as follows 

Xl [n+l]= p iVl [n], < 01 < p? max = Pi/P R ,i, (3) 

where Prj is the received power at the node i and choice of the scaling factor satisfies the power constraint 

In layered networks, all copies of a source signal and a noise symbol introduced at a node traveling along different paths 
arrive at the destination with the same respective time delays. Therefore, the outputs of the source-destination channel are free 



of intersymbol interference. This simplifies the relation between input and output of the channel and allows us to omit the 
time-index while denoting the input and output signals. 

Using ([TJ and ([3]), the input-output channel between the source and destination can be written as 
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where Kl, is the set of L-tuples of node indices corresponding to all paths from the source to the destination with path delay 
L. Similarly, Kij : L-i + i, is the set of L — I + l-tuples of node indices corresponding to all paths from the j th relay of ^ th layer 
to the destination with path delay L — l + 1. 

For all the paths between the source s and the destination t, and all the paths between the j th relay of l lh layer to the 
destination t with path delay L — I + 1, we introduce modified channel gains, respectively, as follows 

hs ^ hsi\fii 1 hi 1 i 2 . . . Pi L h-i L i (5) 

(i u ...,i L )£K L 

h>lj = E ■ ■ ■ Phi L hLi L ,t (6) 
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In terms of these modified channel gain^j the source-destination channel in Q can be written as: 

L n t 

y t = h s x s + ^2 E hl i Zl i + Zt ^ 
i=i j=\ 

Problem Formulation: For a given network-wide scaling vector (3 — (fiu)i<i<L.i<i<7ii, me achievable rate for the channel 
in {7]l with i.i.d. Gaussian input is ((3), (4), (6)): 

7(P s ,/3) = (1/2) log (1 +5^), (8) 

where SNR t , the signal-to-noise ratio at the destination t is: 

P h 2 

SNR t — ^ (9) 

The maximum information-rate Ianc(Ps) achievable in a given layered network with i.i.d. Gaussian input is defined as the 
maximum of I(P s ,f3) over all feasible (3, subject to per relay transmit power constraint ([3]). In other words: 

I A Nc(Ps) d = max I(Ps,P) (10) 
P--o<Pl<Pl :max 

It should be noted that Pu, ma x (the maximum value of the scaling factor for i th node in the I th layer) depends on the scaling 
factors for the nodes in the previous I — 1 layers. 

Given the monotonicity of the log(-) function (all logarithms in this work are assumed to be binary, except some cases 
where we explicitly denote the natural logarithm as ln(-)), we have 

/3 op t = argmax I(P Sl f3)— argmax SNR t (11) 

Therefore in the rest of the paper, we concern ourselves mostly with maximizing the received SNRs. 

In |7), we discussed the computational complexity of exactly solving the problem ( fT0| ) or equivalently the problem ( fTT) . 
Further, we also introduced a key result |7, Lemma 2] that reduces the computational complexity of the problem of computing 
Popt by computing it layer-by-layer as a solution of a cascade of subproblems. However, each of these subproblems itself is 
computationally hard for general network scenarios as it involves maximizing the ratio of posynomials fl4) , p5) , which is 
known to be computationally intractable in general 1 15 1. Therefore in this paper, based on the notion of network simplification 
Q, we introduce an approach to reduce the computational complexity of solving each of these subproblems by selecting 
the best subset of ki relays among the set of ni relays in the I th layer. Thus the proposed scheme leads to an exponential 
reduction in the search space to solve each of these subproblems without significantly compromising on the optimal end-to-end 
performance, in terms of additive and multiplicative gaps between the corresponding performances, as we show below. 

In the following, we discuss the performance of this network simplification based approach to approximate the maximum 
achievable ANC rate first, for the Gaussian N -relay networks and then, for a class of symmetric layered networks. 



'Modified channel gains for even a possibly exponential number of paths as in {5} and {6} can be efficiently computed using line-graphs |l3). Further, the 
number of such modified channel gains scales polynomially in the size of the graph being considered. 



Fig. 2. A diamond network with N relay nodes. 



III. Analog network coding in the Diamond Network: Performance of Network Simplification 

Consider the diamond network of Figure [2] We can consider diamond network as a layered network with only one layer of 
relay nodes. Then using |5]), (|6j, and |9]), we can compute the SNR at the destination t for any scaling vector (3 as 

SNR t = ^ k ° tPlhu)2 (12) 

Therefore, using ( fl~0] > the problem of computing the maximum ANC rate for this network can be formulated as 

max SNR t , (13) 

where /3 = (ft, . . . ,ftv) and /3 maa; = (p 1>max p N ,max) with p 2 max = Pi/{h 2 sl P s + a 2 ), i e TV, TV = {1, . . . , TV}. 

Equating the first-order partial derivatives of the objective function with respect to ft,i £ TV, to zero, we get the following 
TV + 1 conditions for local extrema: 

53 hsiPihu = (14) 
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In Q, we show that all points satisfying ( [T4] i lead to the global minimum of the objective function in ( |13) . Further, we 
show that no real solution of the simultaneous system of equations in ( fl"5j ) exists. In other words, no solution of ( fT3| > exists 
where all relay nodes transmit strictly below their respective transmit power constraints. Therefore the global maximum of 
the objective function occurs at one of TV hyperplanes (of dimension TV — 1) defined by f3 u — /3 u>max ,u 6 TV. Denote the 
hyperplane at which the objective SNR t attains its maximum as f3 u * = (3 U * ,max- The details of how the hyperplane f3 u * is 
identified among TV possible candidate hyperplanes are provided in |8|. 

Solving the system of simultaneous equations in ( fl5] l on the (TV — 1) -dimensional hyperplane defined by f3 u * = P u * :max 
results in 

Pi =J—t - R r ,«eTV\{uj (16) 

"it "su* Pu* ,max" , u*t 

Note that {fif } is the set of scaling-factors for the nodes in the set of relay nodes TV\ {u*} that maximizes SNR t on the 
hyperplane fi u * = P u \max- 

Denote the scaling vector (3 at which SNR t attains its maximum on the hyperplane (3 U + — (3 u *,max as P u = {Pi , ■ ■ ■ , Pn )> 
where pf is given by ( fTo} when i ^ u* and pf — P u * lTna x when i = u*. Computing SNR t ((3 u ) by substituting (3 U in 
(TT~2l results in 



SNR t (f3 u ) = (P s /a 2 ) 



(17) 
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Note that SNR t ((3 u ) in ( fT7| i is computed when /3f in ( fT6] > is not constrained by its corresponding upper-bound f3 i7nax . 
Therefore, SNR t ((3 u ) may not be achievable in the actual system. The optimal achievable destination SNR with TV relays 
SNR° P ^ is obtained by evaluating ( fL?) at (3 op t = (/9^ p * , . . . , /3^f ), where each /3° pt satisfies its corresponding upper bound. 
Theorem 1 and Lemma 2 in |8j provide the exact characterization of (3 op t and detail the procedure of obtaining it from (3 U . 
However, SNR t (f3 u ) provides an upper bound on SNR° P ^, that is: 

SNR°% < SNR t (j3 u *) (19) 

However, this upper bound on SNR° P ^ can be tight or arbitrarily loose depending on the system parameters. In particular, as 
the source power P s increases, this upper-bound becomes looser; on the other hand, as the source power decreases, this upper 
bound becomes tighter. 

In fTT) it is observed that at high rate (or high SNR in the present setting) we should be able to approximate the maximum 
achievable rate within an additive gap and at low rate (or low SNR) we should be able to approximate the maximum achievable 
rate within a multiplicative gap. In this paper, we are particularly concerned with the variation of the destination SNR with 
the source power. Therefore for the ease of presentation, we compute the additive gap for large P s (high-SNR regime) and 
the multiplicative gap for small P s (low-SNR regime). 

Let Rn — (1/2) log(l + SNR° P ^) denote the optimal end-to-end ANC rate achieved by using all N relays in the diamond 
network. Similarly, let Rk — (1/2) log(l + SNR° P £) denote the optimal end-to-end ANC rate achieved by using any k relays 
out of available N relays in the diamond network. In the following, we compute an upper bounds on additive gap (i?jv — Rk) 
and multiplicative gap (R^/Rk) for large and small P s , respectively. 

A. An Upper Bound on the Additive Gap 

It follows from the definition of /3j, the scaling-factor for the i th node, that (3f ~ ~. Then, ( fTo*| ) implies that for a 
sufficiently large P s , say Pj, every /3" violates its corresponding upper-bound and in an optimal scaling-vector /3 opt , we have 
P° pt = Pi^max as P s > Pj. Therefore, with N relays, we have for the optimal SNR at the destination: 

\2 
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where h si . max = max, 6A f h si , h itimin = min ieJv - h it , h it , max = max i£Af h it , P™™ ax = min i6A ^ B i>max , and P™£ x ax = 

max^gj^ /3i t max- 

Consider SNR op ^, the optimal achievable SNR at the destination when any subset JC of k relays is chosen from the set Af 
of all available relays. In other words, SNR op ^ is the SNR expression in ( p"2| ) evaluated at (3 opt corresponding to subset JC. 
Then, we have: 

P s k 2 (h s i. m i n P™™ ax hit, m in)' 2 

— ^2 1 i htamax u ^2 ' ' ' 

1 ' fi \Pi,max n it,max) 
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where h si ^ min = min, eA ^ h s 

Therefore, from ( |20] > and ( |2"T] i, we have 



SNU^x f N\/ ' (hsi,maxf^maxhit,max)' 1 1 + ^iPi^max^it^ax) 2 

SNR op l -\kj (h sitmm (3™ n ax h zt ,rnin) 2 1 + N (fi 
W f N V (h 2 sl . max h tt , max ^ 2 



< 1 — 1 a lt "1 = 7^ , . )2 > W 

rarna% \2 ^2 l-\-k(fi™ ax hit max) 2 

where (a) follows from ^'' m °p < ^'• max and 1+N ^^iT y < 1- Note that parameter a\ characterizes the asymmetry 

in the network and in general, a\ > 1. 

Using d22|, we compute the upper-bound on additive gap, Pjv — Rk, with two approaches: direct and indirect. 



Direct Approach: We compute Rjy — Rk as follows: 



Rn - Rk = x log 



1 1 + snrX 



2 ~ 6 1 + SNRf* 



i^_i + (%) 2 ai SNn$ 

2 
1 



< - log r 

" 2 6 i + SNR° t P l 



< ~\og{{N/kf ai ) 
1 



2 In 2 

Indirect Approach: We compute an upper-bound on Rn — by first expressing it as 

Rn — Rk = (Rn — Rn-i) + ■■■ + (Rk+i — Rk) 

and then summing up the upper-bound on each of these partial gaps. Consider (Rk+i — Rk)'- 

1 1 + (1 + 1/kfaxSNR^ 
Rk +1 -R k <- 2 log 1 - snr7 , 



(21n(AT/fc) +lnai) (23) 
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This leads to 

i N ^ o 

R N -R K <—C£- + (N-k)ln ai ) 

i—k 

<^[ 2 ( H N-i-Hk-i) + (N-k)lii ai ], (25) 

where H n is the Harmonic number, lim n _ J . 00 H n ~ Inn + 7 and 7 is Euler-Mascheroni constant fl6) . 

Combining the upper-bounds on iijv — R k computed in |23) and ( |25] l, we have the following upper-bound on the additive 
gap for sufficiently large P s : 

R N -Rk< ^= min{2 ln(iV/fc) + In a u 2(H N . 1 — H k -i) + (N — k) In ttl } (26) 
2 in 2 

B. An Upper Bound on the Multiplicative Gap 

Equation ( fT6| ) implies that for a sufficiently small P s , say P s s , each /3" is within its corresponding upper-bound. Therefore, 
the upper bound on SNR° P N in ([19) is tight. 

Let SNR™g x denote the destination SNR when each relay in the subset JC transmits at the upper-bound of its transmit 
power constraint, i.e. Pi — Pi.max, i € JC. However, all nodes transmitting at the upper-bound of their transmit power constraint 
results in suboptimal performance, as we showed in J6), ||7). Therefore, we have 

SNR° Pt > SNR" laX — — ^ ig/C hsiPi,maxh.it)' 2 

P k 2 (h ■ B min h-4 ■ ) 2 

i + k{(i™* ax h it , max y 



Therefore, from ( |18) , ( fT9) i, and ((27]), we have 
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where parameter a 2 characterizes the asymmetry in the network and in general, ct^ > 1. 
Using ( |28] l, we compute the upper-bound on multiplicative gap, Rm/Rr, as follows: 

R N log(l + SNR opt N ) 



Rk log(l + SNR°f) 

< log(l + f Q2 (l+ 7 /fc)^iV^) 
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- ¥ a2 ( 1 + fc) - T a2(1 + 7) ' (29) 

where (a) follow when SNR op ^ — > and (b) follows from noting that ( 1+z ) io g (i+ z ) — 1 f° r a ^ z >§. 

The results of the previous two subsections are collected in the following lemma that is our first main result. 

Lemma 1: For the Gaussian TV-relay diamond network, the additive and multiplicative gaps between the optimal performance 
of analog network coding obtained with and without network simplification are bounded from above, respectively, as 

R N ~Rk< J r^min{21n(JV/fc)+lnai,2(ff JV _i-fr fc _i) + (iV-fc)liitti} 
I in I 

Rn <r N fij.'y 
R^^T a \ l+ k 

where ot\ and ct2 are two parameters characterizing the asymmetry in the network. 

Example 1: Consider the Gaussian TV-relay diamond network in the symmetric configuration, in which the channel gains 
to the relays are equal (h s i = h^i € AO and the channel gains from the relays are equal (ha — g,i G AO and in general 
h 7^ g. Also assume that the transmit power constraint on each relay is the same, i.e. E\xf\ < P. Note that in this setting 
Pi max = ft 2 = P/{h 2 P s + c 2 ) and 7 = l/(/3g) 2 . Also, the parameters ot\ and a 2 that characterize the asymmetry in the 
network are ay = «2 = 1- Therefore in this setting, using |26| and <j29j we obtain the following upper bounds on additive 
and multiplicative gaps between Rn and Rk when x > 1: 

Rn - Rk < ~(H N -i - fffe-i) ~ log ^ 
In 2 k 

Rn „. N, 

-^< T (l + 7 /*)< T (l + 7 ), 

where (a) holds asymptotically in and fc. It is interesting to note that these bounds exactly coincide with the similar bounds 
obtained from the first principles for this system configuration, as an interested reader can verify. ■ 

IV. Performance of Network Simplification in General Layered Networks 

In this section we analyze the performance of network simplification in a class of symmetric layered networks with one or 
more layers of relays between the source and destination and each relay performing analog network coding on its input signal. 
The particular class of symmetric networks we consider here are defined such that channel gains along all links between the 
nodes in two adjacent layers are equal. This is a generalization of the symmetric diamond network configuration considered 
in Example 1 above. We introduced such networks in [7] and called them "Equal Channel Gains between Adjacent Layers 
(ECGAL)" networks. Figure [3] provides an illustration of such networks. 

For the ease of presentation, consider ECGAL networks where each layer of relay nodes has N relays and all relay nodes 
have the same transmit power constraint EX 2 < P. 

Using the following result we establish in ||7j, we can obtain the optimal solution (3^ pt of problem ( fTl) for an ECGAL 
network with N nodes in each relay layer. 

Lemma 2 (Lemma 2 in Consider a layered relay network of L + 2 layers, with the source s in layer '0', the destination 
t in layer 'L + V, and L layers of relay nodes between them. The I th layer contains n; nodes, Uq = = L A network- wide 
scaling vector f3 op t = (/3° pt , ■ ■ ■ ,/3i, pt ) that solves (JTTJ for this network, can be computed recursively for 1 < I < L as 

ni+i 

I3° pt = argmax TT (1 + SNR l+u {(3{ p \ ft*, fl)), 

where (3° pt is the subvector of optimal scaling factors for the nodes in the / th layer, (3° pt — (f$° p , . . . and constraints 

A 2 < Ptmax ^ component-wise /3 2 < fil max . ■ 




Fig. 3. An ECGAL network of L + 2 layers, with the source s in layer '0', the destination t in layer 'L + 1', and L layers consisting of two relay nodes 
each between them. The channel gains along all links between two adjacent layers are equal. 

The solution (3^ t has the property that under it all relays in a layer use the same scaling factor and it is equal to the 
maximum value of the scaling factor for the nodes in the layer, i.e. for each layer i, 1 < I < L, and every node i in the I th 
layer, 1 < i < N, (/^) 2 = (/^) 2 = where 
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Therefore, (3^ pt can be written as j3^ pt = 
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* A "" K " "" "* A ^ ' <A . . ,(3^). We use superscript JV to emphasize that the optimal scaling factors 

in p0| ) are computed for nodes in each layer. 

The optimal SNR at the destination with N relays in each layer of a layered network with L relay layers is given as follows 

L \ 2 

1=1 
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where /3j is given as in <[30j>. 

Now consider the network simplification scenario where only k out of N available relays in each layer are used. Using Lemma 
2 in |7j, we can solve problem ( fTTj ) for this simplified ECGAL network and obtain the optimal solution 0% pt — (/3±, . . . , /?£), 
where fi\ is the optimum scaling factor for all nodes in layer I and is given as follows: 
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The corresponding optimal SNR at the destination is given as: 



SNR°X = 



P, 



L 

n 

1=1 



t,k 



(33) 



Using ( (31} and ( |33] ), we can compute the ratio of the optimal destination SNRs obtained when all N relays in every layer 
are used and when only k out of N relays in every layer are used. After a few steps of algebraic manipulation, we get 



SNR opt N 
SNR° P * 



2L-X 



*=i ( n (kPfhj)) 2 (nikpfh,))' 



(34) 



Ett^ 



i=i(nwf(.j)) ! (nww) : 



Let Rn = (1/2) log(l + SNR° P ^) denote the optimal end-to-end ANC rate achieved by using all N relays every layer of 
an ECGAL network. Similarly, let Rk — (l/2)log(l + SNR° P £) denote the optimal end-to-end ANC rate achieved by using 



any k relays out of available N relays in every layer of an ECGAL network. Arguing as in Section III in the following we 
compute the additive gap Rjy — Rk for large P s and the multiplicative gap R^/Rk for small P s . 

Before we discuss the upper-bounds on the additive and multiplicative gaps for ECGAL networks with arbitrary number of 
relay layers, consider the following example where we compute such bounds for an ECGAL network with two layers of relay 
nodes (L = 2) for any N and k. 

Example 2: Consider an ECGAL network with two layers of relay nodes between the source and the destination, L = 2. 
Using ([34]), we have for this network 



SNR°% 
SNR° t p * 



1 



1 + 



1 



1 



k 2 hjpf k^h\h 2 2 l3l{^f 



+ 



N 2 h\j3l N^hlhlPlifi^) 2 
Substituting for and fi\ from ( |3T| i and ( [33] l in the above expression results in 



(for P s 



oo 



(for P s -> 0) 



SNR°% 
SNR°$ 

SNR°% 
SNR° P * ~~ 



4 1 



k) 1 



K 

khj 



k 2 h 2 P 



3 1 



hi 
Nh : 



k 2 P 



■■ ^ N 2 h 2 P 



1 



(35) 



(36) 



Nh 2 h 2 P) 



N 2 p\iq 

Following the same sequence of steps as in the direct approach to compute the upper bound on the additive and multiplicative 



gaps for the diamond network in subsections III-A and III-B respectively, using (|35fl and (pojl, we get the following upper 



bounds on the corresponding additive and multiplicative gaps for the two layer ECGAL network: 

1 



Rn — Rk < 



1 ^ N 

2lrT2 41n T 



In- 



hi 

khj 



k 2 h 2 P 



1 



h 2 
Nh 2 



N 2 h 2 P 



R 



N 



<(- 

Rk ~ \ k 



3 1 



k 2 P 



1 -4- -^—( 1 -I- 1 



for arbitrary N and k relays in each layer. ■ 
However, for ECGAL networks with arbitrary number of relay layers, it is analytically hard to compute such upper-bounds 
on the additive and multiplicative gaps between the optimal end-to-end performances with and without network simplification 
for any N and k. Therefore, in the following we attempt to analyze the scaling behavior of such upper bounds with large N 
and k. 

SNP° pt . i—i 

Consider the asymptotic behavior of - 'ff with P s and large N, and k. Using ( f30] > and ( |32| ), we have for large N and k: 

bJyH tk 



(for P s -t oo) 
(for P s ->■ 0) 



ri(A 

»=i 



N\2 



P/(h 2 s P s ) 



i=l 
I 



Then, from ( [34] > and ( [37] i we obtain for large N and k, and P s — > oo: 

2L r 



=i 



KKL 



KKL 



SNR° t % 



, 1 1 

1 + a 

k N 



SNR% 

and from ( f34| > and ( |38| ) we obtain for large N and fc, and P s — > 

1 



l-x 1 

where a = h\ ^ ^ 



SNR° t % 
SNRfi 



2L-1 



1 + 6 



1 

iV 2 



where b = 



P \hl 



hi 



(37) 
(38) 

(39) 

(40) 



With these results we are now ready to compute the upper bounds on additive and multiplicative gaps between the optimal 
performance of analog network coding with and without network simplification. 



First, we consider the upper bound on the additive gap. Following the same sequence of steps as in the direct approach to 



compute the upper bound on the additive gap for the diamond network in subsection III-A using ( |39[ l we obtain the following 
upper bound on the additive gap Rn — Rk for P s — > oo and asymptotically large N and k satisfying N,k — o(P s ): 



R 



N 



R k < 



1 



N 

2Lln — 
k 



In 



1 



a ~ - — 



(41) 



2 In 2 

Next, we consider the upper bound on the multiplicative gap. Following the same sequence of steps as in subsection |III-B| to 
compute the upper bound on the multiplicative gap for the diamond network, using ( |40| > we obtain the following upper bound 
on the multiplicative gap Rn /Rk for P s — > and asymptotically large N and k satisfying N,k = o(l/P s ): 

2L-1 r 



Rn 

Rk 



< 



1 



1 

A 3 



1 

AT2 



(42) 



These results on the asymptotic behavior of the additive and multiplicative gaps are collected in the following lemma that 
is the second main result of this paper. 

Lemma 3: For ECGAL network, the asymptotic additive and multiplicative gaps between the optimal performance of analog 
network coding obtained with and without network simplification are bounded from above as 



Rn — Rk < 



1 



2 In 2 
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where a = h\ J] 



Rn 

L ~l 1/hf and b = 
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N 



ln(l 



1 

k 2 



l 

n 2 



< 



2L-1 



-{l/h\ + are constants depending on the system parameters. 



Remark 1: These bounds imply that each relay layer increases the additive gap by 2\og(N/k) and the multiplicative gap 
by a factor of (N/k) 2 . Thus the per layer contribution to the gaps is in agreement with the corresponding results obtained 



for the diamond network (with a single relay layer) in Section III For instance, for the network in Example 1 in Section III 
using the results of this section we obtain 



R N -Rk< \og{N/k) 



R_N_ < N_ 
Rk ~ k 



1 



Pa 1 



with L = 1, a = and hi = g. These bounds coincide exactly with the corresponding upper-bounds obtained in Example 1. 
This indicates that the bounds in ( |4"Tj ) and ( |42| ) are actually tight. 

Remark 2: The asymptotic expressions of Il^iC/^fO 2 an d Y[i=x(Pi) 2 m 'HZr an< ^ *HD com P utec l f° r large N and k result 
in smaller values of the ratio ' '™ t computed in d39b and (HO) compared to the corresponding values computed using the 

expressions for and f3^ in ( f30| > and < [32] >, respectively, for any N and k. This implies that the upper bounds on the additive 
(RN — Rk) and multiplicative (Rn / Rk) gaps in ( pT| and |42) computed for asymptotically large N and k are smaller compared 
to the corresponding bounds Rn — Rk and RM/Rk obtained for any N and k. However, it can be proved that Rn — Rk 
and Rn / Rk approximate Rn — Rk and RN/Rk, respectively, within constants, depending on various system parameters, i.e. 
Rn — Rk < Rn — Rk + c and jg- < ^-(1 + d), where c and d are two constants dependent on system parameters. For 



instance, for the network in Example 2, we can show that Rn — Rk < Rn 

V. Conclusion and Future Work 



i? fc + ^and 



Computing the maximum end-to-end rate achievable with analog network coding in general layered relay networks is an 
important but computationally intractable problem. Previously, this problem is addressed under certain assumptions on input 
signal factors and received SNRs at the relay nodes. However, such assumptions lead to suboptimal end-to-end performance 
in general, and in certain communication scenarios it is desirable to compute the optimal rates without such assumptions. 

We introduce an approach based on the notion of network simplification to approximate the optimal ANC rate within small 
additive and multiplicative gaps in the Gaussian A^-relay diamond network and a class of symmetric layered networks while 
simultaneously reducing the computational complexity of solving this problem. One of our main contributions is a result that 
states that in a layered network, each layer increases the additive gap between the optimal performance of analog network 
coding without and with network simplification by log(iV /k) 2 and similarly, the multiplicative gap by (N/k) 2 . To the best of 
our knowledge, this work provides the first characterization of the performance of network simplification in general layered 
amplify-and-forward networks. Also, our results suggest a new approach to approximate the optimal ANC rates while allowing 
for the computation of the additive and multiplicative gaps simultaneously. In future, we plan to extend this work to more 
general networks and to refine the upper bounds on the additive and multiplicative gaps computed in this paper. 
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